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Abstract. We investigate the non-equilibrium thermal quantum discord and 
entanglement of a three-spin chain whose two end spins are respectively coupled to two 
thermal reservoirs at different temperatures. In the three-spin chain, besides the XX- 
type nearest-neighbor two-spin interaction, a multi-spin interaction is also considered 
and a homogenous magnetic field is applied to each spin. We show that the extreme 
steady-state quantum discord and entanglement of the two end spins can always be 
created by holding both a large magnetic field and a strong multi-spin interaction. 
The results are explained by the thermal excitation depression due to switching a 
large energy gap between the ground state and the first excited state. The present 
investigation may provide a useful approach to control coupling between a quantum 
system and its environment. 
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1. Introduction 

Quantum entanglement was once considered as a unique resource that can be used 
in quantum information processing pQ. However, recent researches have shown that 
besides entanglement a composite quantum system may have other kinds of nonclassical 
correlations which can appear even in separable states [21 El HJ El E]- In order to 
quantitatively describe such quantum correlations in a composite quantum system, 
many different measures have been proposed P, El [9] . Among them, quantum discord 
(QD), firstly introduced by Ollivier and Zurk [6], has received considerable attention 
[TOj EH [121 1131 EH EEl EE] • Using an optical architecture, Lanyon et al. [I] experimentally 
showed that even fully separable states with quantum discord can be used to construct 
quantum computer. 

In all real situations, quantum systems can not be completely isolated from their 
environments. Coupling of a quantum system to its surrounding unavoidablely results 
in the destruction of quantum correlations of the system. The effect of environments 
to entanglement of bipartite quantum systems have intensively been investigated. It 
has been shown that entanglement undergoes sudden death due to the interaction 
of quantum systems with their reservoirs [TTj . In recent years, the QD dynamics of 
open quantum systems has also attracted much interest in both theory and experiment 
[181II91EQ1EI1E21E31I21E51I261E7]. Werlang et al. p] investigated the dynamics 
of both entanglement and QD in the Markovian environments, and showed that QD is 
more robust against decoherence than entanglement. Recently, an interesting dynamical 
feature of QD, named sudden transition, has been observed [2TI [22] . It means that for 
certain initial states QD undergoes sudden change between a "classical decoherence" 
phase and a "quantum decoherence" phase [2D1 EI] • This sudden transition behavior can 
be explained in a geometrical way and has connection with the property of environment 
[201 12T] . Xu et al. (2E1 EE] experimentally investigated both the Markovian and non- 
Markovian dynamics of classical and quantum correlations and observed the sudden 
transition behavior of QD. 

Apart from the situation in which a quantum system is coupled to a single 
environment, it may also be possible that a quantum system is simultaneously in 
contact with two different thermal baths. In semiconductor quantum dots nuclear spins 
and electronic spins consist of a composite quantum system for quantum information 
processing and quantum computing but the coupling manners of nuclear spins and 
electron spins to their surroundings are much different [211 [221 [30]. With the help 
of NMR and quantum optical techniques, one can create two reservoirs at different 
effective temperatures for nuclear spins and electron spins in quantum dots [31], [32] . For 
superconductor qubits, the two-different-thermal-bath coupling situation may directly 
be designed [33]. When interacting with two reservoirs at different temperatures, 
a quantum system may approach a steady state instead of a thermal equilibrium 
state. Thus, in general, the presence of heat /energy /mass currents passing through 
the quantum system in a steady-state may modify the quantum correlations. 
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In recent years, quantum correlations of coupled qubits in contact with two 
different thermal environments, i.e., the non-equilibrium thermal environment model, 
have received some attention [311 EHJ EHJ [371 EE]- Quiroga et al. [36] proposed 
a two-interacting spin- 1/2 system in contact with two heat reservoirs at different 
temperature, and identified a nonequilibrium enhancement-suppression transition 
behavior of entanglement due to the presence of temperature gradient. Employing the 
same model, Sinaysky et al. [37] found that the spin system can converge to steady state 
and studied the dependence of the steady-state concurrence on the mean temperature 
and temperature difference of the reservoirs and the energy splits of the spins. Huang 
et al. [38] investigated the nonequilibrium thermal steady-state entanglement in a three 
spin- 1/2 XX chain in contact with two heat reservoirs at different temperature and 
found that the temperature difference of the heat baths benefits the entanglement in 
the nonsymmetric coupling case. Spin chain models in contact with two reservoirs at 
two different temperatures have also been employed for studying heat current transfer 
[391 HOj HI] . Yan et al. [39] considered an interacting spin- 1/2 chain connected to two 
phonon baths held at different temperatures and showed that heat transport through 
the spins systems can be controlled by applying an inhomogeneous magnetic field due 
to switching an energy gap. 

Stimulated by the pervious investigations, we here study how to control the steady- 
state QD and entanglement of spin systems. We consider a three-spin- 1/2 chain in which 
besides the XX-type nearest-neighbor two-spin interaction a three-spin interaction is 
included and an external magnetic field is homogeneously applied to each spin, and 
meanwhile the two end spins are coupled to two thermal environments at different 
temperatures. We show that the coupling of the spin system to the thermal reservoirs 
can be controlled and the thermal excitation can be greatly depressed by the three- 
spin interaction and the magnetic field. As a result, the extreme steady-state QD and 
entanglement in the two end spins can be created. 

This paper is organized as follows. In Sec. II, we describe the model and introduce 
the calculation method. In Sec. Ill, definitions on quantum discord and concurrence 
are briefly reviewed. In Sec. IV, numerical results, discussion and physical explanations 
are presented. Finally, a summary is given in Sec. V. 

2. Model And Master Equation 

The model under investigation is described in Fig. 1. We consider a three-spin- 1/2 
chain which Hamiltonian reads 

Hs = JY, («• , + + h E < + k K^s + ofoM) , (1) 

i=l i=l 

where erf (a = x, y, z) are the Pauli matrices for the zth spin, J is the coupling constant 
between the nearest-neighbor spins, and h is the external magnetic field strength, 
homogeneously applied to each spin. Besides the two-spin interaction, the three-spin 
interaction [121 E3 HI] is also included, which strength is denoted by k. 
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Figure 1. A schematic representation of a three-spin chain coupled to two thermal 
baths at different temperatures, Ti and T 3 . 



As shown in Fig. 1, two end spins 1 and 3 are in contact with two phonon 
baths at different temperatures, T\ and T3, respectively. In the interaction picture, 
the Hamiltonian describing the interaction between the jth spin and its phonon bath is 
given by 



H 



SBj 



J.) 



a*®B j ,{j = 1,3), (2) 



where b\j{b n j) is the creation (annihilation) operator for the nth mode of thermal bath 
j, and gj is the coupling constant between the jth spin and the nth bath mode. 
Let us first consider the eigenvalue problem of the Hamiltonian (1) 

H s \<j>i) = e l \<j> l ),(! = l,2,3,...,8). (3) 

The spin-up and spin-down states of spin i are represented by state- vectors and 
1 ) ^ , respectively. In the presentation spanned by the uncoupled basis |nin 2 n 3 ) = 
|ni) 1 <S> ^2)2 ® 1^3)3 with rii = 0,1, we can easily work out the eigenstates of Eq. 



as follows 
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1 1 
|0 6 ) = — r sma 2 |110) - cosa 2 1 101) + -^sina 2 |011) , (9) 
v2 v2 

|0 7 ) = -isina 2 |100) + cosa 2 |010) + 4=sina 2 |001), (10) 
\/2 v2 

|0 8 ) = -=sinai|110) -cosai|101) + -=sinai|011), (11) 
v 2 v 2 

with the corresponding eigenvalues £i = — 3/i, e 2 = 3h, £3 = h — 2/c, e 4 = — /i + 2k, £5 = 

—h — k — B, eq = h + k — B,Er = —h — k + B, £ 8 = h + k + B, where B = \/8 + A; 2 , sin cki = 



2V2/J8 + (k - B)\ cosai = (k — B) j J8 + (k — B) 2 , sina 2 = 2V2/V8 + (k + B) 2 



and cos a 2 = (k + B) /^8 + (k + B) 2 . 

In the representation spanned by eigenstates f|4|)- (|TT|) . the Hamiltonian of the 
coupled reservoir -spin system can be written as 

H = H s + H SBX +Hsm = Y, e i\<t>i)(<t>i\+ HH A M® B v ( 12 ) 

1=1 5=1,3 w 

where 

A J (co) = (0il^l0i'>l0i>(0/'l- (13) 

In Eq. ffl2l) . the summation XL m ust be done over all possible differences between any two 
eigenenergies of the Hamiltonian (1). In Eq. ( ITBl) . the summation Z) e; - e;; =aj is over a U 
the eigenvalues with a fixed difference u. Obviously, A\u) = A(—cu). Upon substitution 
of eigenstates (l4l-(fTTT) into Eq. ( fl3l) . we find the following nonzero transition operators 

A (wi) = -^=(sinai |0 2 ) (4>s\ ~ cosa 2 \<j> 6 ) (0 4 | 

- cosa 2 |0 3 ) (07 1 + sinai |0 5 ) (0i|), (14) 

A (w 2 ) = -^(sina 2 |0 2 ) (06 1 - coscti |0 8 ) (0 4 | 

- cosai |0 3 ) (0 5 | + sina 2 |0 7 ) (0^), (15) 

A (w 3 ) = 4= (1 02) (03 1 - sin (a_) |0 6 ) (0 5 | 

+ sin(a_)|0 8 >(0 7 |-|0 4 >(0 1 |), (16) 

^3 ( w i) = ~j^( sinai 1^) (08 1 + cosa 2 

+ cosa 2 |0 3 ) (0 7 | + sinai |0 5 ) (0i|), (17) 

^3 (^2) = -^(sina 2 |0 2 ) (0 6 1 + cosai |0 8 ) (0 4 | 

+ cosai |0 3 ) (05 1 + sina 2 |0 7 ) (0i|), (18) 

A M = ^t=(|02> (03 1 - sin (a+) |0 6 ) (0 5 | 

-sin(a + )|0 8 )(0 7 |-|0 4 )(0 1 |), (19) 
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where W\ = 2h — k — B, uj 2 = 2h — k + B, u 3 = 2(h + k), a + = a,\ + a 2 and a_ = a\ — a 2 . 

By means of the general reservoir theory within the Born-Markov and rotating 
wave approximations [321 EES HZ] , one can obtain the equation of motion for the reduced 
density matrix of the spin chain 

^^-i[H s ,p)+L 1 (p)+L 3 (p), (20) 

where Lj(p) (j — 1, 3) is the dissipative term due to the coupling of spin j to its thermal 
bath and is given by 

Hp) = E 7iK)(l + n^)) (2Aj (w M ) pA) (w M ) - {p, A] (^) ^ 

w,,>0 

+ E 7j(w>j(w„) (2A} {uj pAj - {p, Aj (^) At (^)})(21) 

o; Al >0 

In deriving out the master equation ( J20l) . we have assumed that the jth bath 
is always in a thermal equilibrium state at temperature Tj. In Eq. f)2ip . nj{u^) = 
l/(exp(/3jU ll ) — 1) with 0j — 1/ (Tj) is the mean thermal photon number of the jth bath 
at frequency (taking the Boltzmann constant ks — 1), and 7j(w M ) is defined through 
the integral vrX)„ Ifl 1 ]™^ 2 (l+ < b^b n j = J °° 7^(0;^) (1 + n J (a; /1 ))da; At . Here, the Lamb 
shift has been omitted. 



3. Quantum Discord and Concurrence 

In this section, for convenience of discussions in the next section, we give a brief 
review on quantum discord (QD) and concurrence. QD is defined as the discrepancy 
between quantum extensions of two equivalent expressions for the classical mutual 
information [6J. In classical information theory (CIT), the total correlation between two 
random variables A and B can be described by either the mutual information [TJ HHl H9] 

I C (A:B) = H(A) + H(B)-H{A,B) (22) 

or the equivalent expression 

J C {A:B) = H{A)-H{A\\B), (23) 

where H(X) = —J2xPx^og 2 p x (X = A, B and AB) is the Shannon entropy of the 
variable X with p x being the probability of X assuming the value x, and H(A\\B) = 
-T,a,bPab^og 2 p a \ b = H(A, B) - H(B) (p a \ b = Pab/pb) is the conditional entropy, which 
represents a weighted average of the entropies of A given the value of B. 

In the quantum information theory (QIT) (TJ HH [49] , the total correlation of a 
bipartite system consisting of subsystems A and B in a state described by the density 
matrix pab is defined as 

h (pa-.b) = S(p A ) + S(p B ) - S(p AB ), (24) 

which is the straightforward extension of ( |22i) . Here, S(pa(b)) — — Tt(pa(b) l°g2 Pa{b)) 
is the von Neumann entropy of the subsystem A(B), while S(pab) = ~Tt(pab log 2 Pab) 
is the entropy of the composite system AB. 
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The extension of (122 p to the quantum realm is no longer straightforward since the 
value of H(A\\B) is measurement dependence, and quantum measurement may fully 
destroy a quantum state. The counterpart of ( 1231 in QIT may be defined as 

J q (pa-.b) = S(p A ) - S {nf} {pa\b) , (25) 

where jn^j are a set of projectors performed locally on subsystem B, and 
S{ n Bj \Pa\b) = 'EjqjS{p :1 A ) with p\ = Tr B (nfp AB Tlf) /qj and the probability qj = 
Tr A J B {Ufp AB Uf). The project operator Uf = \6j) (6j\ with |^) = cos0 |O) + e^sin0 |1) 
and \0 2 ) = -cos6 |1) + e-^sin0|O) (0 < 6 < 2tt, < <f) < 2tt). From $25]), it is clear 
that different choices of jn^j may lead to different values of J q (pa-.b)- The minimum 
difference between I q (pa-.b) and Jq {.Pa-.b), called quantum discord (QD) [6], is used to 
describe the quantum correlation of a bipartite quantum system 

D (Pa-.b) = pin [I q (pa-.b) ~ Jq (Pa-.b)} (26) 
or equivalent ly 

D (Pa-.b) = Iq (pa-.b) - max [J q (pa-.b)} ■ (27) 

{ n f} 

From ( 1241) and ( 1271) . the classical correlation contained in a quantum system is 
defined as [2] 



(28) 



C (pab) = I q (pa-.b) - D (Pa-.b) = max S (p A ) - Sr nB \ (pa\b) 

{ n f } L 1 J 1 V 1 

In our investigation, entanglement is qualified by the Wootters concurrence |50j . 

For given density matrix pab of a bipartite system AB , the concurrence is defined as 

C = max jo, y/Xi — y/Xz — — v/A^}, where A« (z = 1,2,3,4) are the eigenvalues 

of the matrix R = p ab {p\ ® ° V b) P*ab ( a A ® ° V b) > arranged in decreasing order of 

magnitude, is the complex conjugate of p^s and B are the Pauli matrices for 

systems A and B. The concurrence attains its maximum value 1 for maximally entangled 

states and for separable states. 



4. Results and Discussion 



The master equation (1201 can be easily solved by the fourth-order Runge-Kutta 
method in the representation spanned by the eigenstates of Hs- We take the evolution 
time long enough such that the final density matrix reaches steady state p s t- Then, 
according the definitions on QD and concurrence given in the preceding section, we can 
investigate the influence of the bath temperature, multi-spin interaction and external 
magnetic field on the QD and concurrence of the spin chain. In the calculation, we 
set the coupling constant ,1—1. It means that all the interaction constants in the 
Hamiltonian are rescaled by the XX spin chain coupling strength. We also assume that 
the decay rate is spectrum independent, i.e. 7(0;) = 7. 
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Figure 2. Steady-state QD (green lines) and concurrence (red lines) as a function of 
the field h with various values of the temperature difference. The other parameters 
are chosen to be 7 = 0.01, Jm = 1-8 and k — 2. The figures (a) and (b) are for the 
spin pairs 13 and 23, respectively. 




Figure 3. Steady-state QD (green lines) and concurrence (red lines) as a function of 
the field h with various values of the mean temperature. The other parameters are 
chosen to be 7 = 0.01, AT = 0.5 and k = 2. The figures (a) and (b) are for the spin 
pairs 13 and 23, respectively. 



In Figs. 2 and 3, the steady-state QD and concurrence of spin pairs 13 and 23 are 
shown as a function of the magnetic field for various values of the temperature difference 
AT = T\ — T 3 and of the mean temperature Tm = (T\ + T 3 )/2. In these figures, we 
see that both the QD and concurrence first increase with increasing of the field, get 
maximal values and then decay to zero. As either the temperature difference or the mean 
temperature increases, in general, the QD and concurrence are diminished. Comparing 
Fig. 2 with Fig. 3, we notice that the mean temperature affects more strongly the 
concurrence than the temperature difference. The sudden death of concurrence as 
shown in Fig. 3 indicates that the QD is more robust against the mean temperature 
than the concurrence. From these figures, we come to the conclusion that the QD and 
concurrence can be enhanced by switching on the properly large magnetic field if both 



the mean temperature and temperature difference are not large. 
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Figure 4. Steady-state QD (green lines) and concurrence (red lines) as a function of 
the three-spin interaction strength k with various values of the mean temperature and 
of the field strength. The figures (a) and (c), and (b) and (d) are for spin pairs 13 and 
23, respectively. The other parameters are chosen to be 7 = 0.01 and AT = 0.8. The 
symbols shown in the inset of Fig. 4(b) are applicable to all the curves of Fig. 4. 

Figures 4 and 5 show the steady-state QD and concurrence of spin pairs 13 and 23 
as a function of the three-spin interaction strength k for various values of the magnetic 
field strength. In these figures, we see that the QD and concurrence for spin pair 13 
first increase with the three-spin interaction and then get a plateau. It is very interest 
that the plateau can be raised to the maximum level by increasing the magnetic field 
strength. As shown in Figs. 4 and 5, this feature can be maintained even if the mean 
temperature is high and the temperature difference is large. As for the spin pair 23, 
figures 4 and 5 show that its QD and concurrence first increase with the three-spin 
interaction, get peaks and then decay to zero. In Fig. 6, the steady-state QD and 
concurrence of spin pairs 13 and 23 as a function of the magnetic field strength h with 
various values of the three-spin interaction strength k. In Fig. 6(a), we observe that the 
QD and concurrence of spin pair 13 can get the maximum level plateau by increasing 
the field if the three-spin interaction is enough strong. The maximum plateau width 
is enlarged as the interaction strength increases. Thus, we can maintain the extreme 
QD and concurrence of the spin pair 13 by holding the strong interaction and magnetic 
field. 

In order to find out the physical reasons for the observed phenomena, we first 
analyze the eigenvalues ei (I = 1,...,8) of the Hamiltonian (1) as a function of the 
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Figure 5. Steady- state QD (green lines) and concurrence (red lines) as a function of 
the three-spin interaction strength k with various values of the temperature difference 
and of the field strength. The figures (a) and (c), and (b) and (d) are for spin pairs 
13 and 23, respectively. The other parameters are chosen to be 7 = 0.01 and Tm = 2. 
The symbols shown in the inset of Fig. 5(b) are applicable to all the curves of Fig. 5. 
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Figure 6. Steady- state QD (green lines) and concurrence (red lines) as a function of 
the magnetic field strength h. The figures (a) and (b) are for spin pairs 13 and 23, 
respectively. The parameters are chosen to be 7 = 0.01, Tm — 1.2 and AT = 0.8. The 
symbols shown in the inset of Fig. 6(b) are applicable to the curves of Fig. 6 (a). 
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Figure 7. Eigenenergy Ei of the Hamiltonian (1) as a function of the three-spin 
interaction strength k with various values of the held strength. The symbols shown in 
the subset of Fig. 7(a) are applicable to all the curves in the figures. 




Figure 8. The energy difference between £ 3 and £5 as a function of the field with a 
fixed value of the interaction, k = 10. 
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interaction and magnetic field strengths. When the field is weak(/i << k), the eigenstate 
|05) is the ground state of the spin chain since the eigenenergy £ 5 = —h — k — + k 2 is 
the smallest one as shown in Fig. 7(a). It is noted that E\ = —3h. Thus, the eigenstate 
\4>i) becomes the ground state when h > k. The two states have the energy crossing 
around the point h = k as shown in Figs. 7 (b)-7(d). Since then, the eigenstate 1 0s) 
becomes the ground state. As the interaction strength k further increases, the state \<ps) 
with the eigenenergy e 3 = —h — 2k crosses with the state \<pi) and becomes the first 
excited state of the spin chain. In Fig. 8, the energy difference between the eigenstates 
\(f> 5 ) and |0 3 ) is plotted as a function of the magnetic field. We see that the energy 
splitting linearly increases as the magnetic field increases. In fact, we have £3 — £5 « 2h 
when k is large. 




5 10 15 20 10 20 30 40 
k k 



Figure 9. Eigenstate occupation probabilities of the steady-state density matrix p st 
as a function of the three-spin interaction strength k with various strength of the field. 
The other parameters are chosen to be 7 = 0.01, Tm — 1.2 and AT = 0.8. The symbols 
shown in the subset of Fig. 9(a) are applicable to all the curves of Figs. 9. 

Figure 9 shows the eigenstate occupation probabilities which are defined as Pi = 
tr(\<pi){4>i\pst) = ((j>i\Pst\4>i) as a function of the three-spin interaction. We see that is 
the most populated state when k « h, \<j>i) and |0 5 ) cross and take the same probability 
around the point k = h when h is large, and then becomes the most populated 
state and takes over all the occupation probability, as shown in Figs. 9(b)-9(d). 

Therefore, when the magnetic field and three-spin interaction are strong, the 
possible thermal excited transition is one from to [fa), which is induced by the 
thermal resources interacting with the spins 1 and 3, respectively. However, E3—E5 ~ 2h, 
as discussed above. Thus, the thermal excitation can be mostly depressed if h ^> T\, T 3 . 
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Figure 10. The QD and concurrence of \4>s) as a function of the three-spin 
interaction with various values of the field. The other parameters are chosen to be 
7 = 0.01, Tm = 1.2 and AT = 0.8. The curves (a) and (b) correspond to correlations 
of spin pairs 13 and 23, respectively. The symbols shown in the subset of Fig. 10(b) 
are also applicable to all the curves of Fig. 10(a). 

In this case, the spin chain is nearly decoupled from the thermal resources. Thus, the QD 
and concurrence of the spin system are determined by the most populated eigenstate, 
i.e. the ground state of the spin system. These results mean that the eigenstate |0s) 
makes the most contribution to the QD and concurrence of the spin pairs 13 and 23 
when the magnetic field and three-spin interaction are strong enough. In Fig. 10, the 
QD and concurrence of |0s) are plotted as a function of the three-spin interaction k. 
Comparing Fig. 10 with Fig. 4, we see that the above conclusion is really true. 

5. Summary 

We investigate the quantum discord and concurrence of the three-spin chain which 
ends are coupled to two thermal reservoirs at different temperatures. Besides the XX- 
type nearest-neighbor two-spin interaction a three-spin interaction is also included and a 
homogenous magnetic field is applied to each spin. For fixed temperatures of the thermal 
reservoirs, we find that the extreme steady-state QD and concurrence of the two end 
spins can always be created by raising the magnetic field strength with a strong multi- 
spin interaction. We show that the energy gap between the most populated ground state 
and the first excited state of the spin chain can become much larger than the thermal 
excitation energy when the magnetic field and multi-spin interaction are strong enough. 
In this way, the thermal excitation induced by the thermal reservoirs is nearly depressed 
and the spin chain is decoupled from the thermal environments. As a result, the QD 
and concurrence of the spin chain are totaly determined by the most populated ground 
state of the spin chain. The present results may provide a useful approach to control 
coupling between a quantum system and its environment. 



14 



Acknowledgments 

This work was supported by the National Basic Research Program of China 
(973 Program) No. 2010CB923102, Special Prophase Project on the National Basic 
Research Program of China (Grant No.2011CB311807) and the National Nature Science 
Foundation of China (Grant No. 11074199). 

References 

[1] Neilsen M A and Chuang I L 2000 Quantum Computation and Quantum Information (Cambridge: 

Cambridge University Press) 
[2] Henderson L and Vedral V 2001 Classical, quantum and total correlations J. Phys. A 34 68899 
[3] Vedral V 2003 Classical correlation and entanglement in quantum measurements Phys. Rev. Lett. 

90 050401 

[4] Lanyon B P, Barbieri M, Almeida M P and White A G 2008 Experimental quantum computing 

without entanglement Phys. Rev. Lett. 101 200501 
[5] Knill E and Laflamme R 1992 Power of one bit of quantum information Phys. Rev. Lett. 81 5672 
[6] Ollivier H and Zurek W H 2001 Quantum discord: a measure of the quantumness of correlations 

Phys. Rev. Lett. 88 017901 
[7] Oppenheim J, Horodecki M, Horodecki P and Horodecki R 2002 Thermodynamical approach to 

quantifying quantum correlations Phys. Rev. Lett. 89 180402 
[8] Luo S 2005 Using measurement-induced disturbance to characterize correlations as classical or 

quantum Phys. Rev. A 77 022301 
[9] Modi K, Patarck T, Son W, Vedral V and Williamson M 2010 Unified view of quantum and 

classical correlations Phys. Rev. Lett. 104 080501 
[10] Maziero J, Celeri L C, Serra R M and Vedral V 2009 Classical and quantum correlations under 

decoherence Phys. Rev. A 80 044102 
[11] Dillcnschneider R 2008 Quantum discord and quantum phase transition in spin chains Phys. Rev. 

B 78 224413 

[12] Sarandy M S 2009 Classical correlation and quantum discord in critical systems Phys. Rev. A 80 
022108 

[13] Werlang T and Rigolin G 2010 Thermal and magnetic quantum discord in Heisenberg models 
Phys. Rev. A 81 044101 

[14] Werlang T, Trippe C, Ribeiro GAP and Rigolin G 2010 Quantum correlations in spin chains at 

finite temperatures and quantum phase transitions Phys. Rev. Lett. 105 095702 
[15] Giorda P and Paris M G A 2010 Gaussian quantum discord Phys. Rev. Lett. 105 020503 
[16] Adesso G and Datta A 2010 Quantum versus classical correlations in Gaussian states Phys. Rev. 
Lett. 105 030501 

[17] Yu T, Eberly J H 2004 Finite-time disentanglement via spontaneous emission Phys. Rev. Lett. 93 

140404; Yu T, Eberly J H 2006 Achieving peak brightness in an atom laser ibid. 96 140403; Yu 

T, Eberly J H 2009 Sudden death of entanglement Science 323 598 
[18] Werlang T, Souza S, Fanchini F F and Villas Boas C J 2009 Robustness of quantum discord to 

sudden death Phys. Rev. A 80 024103 
[19] Fanchini F F, Werlang T, Brasil C A, Arruda L G E and Caldeira A O 2010 Non-Markovian 

dynamics of quantum discord Phys. Rev. A 81 052107 
[20] Mazzola L, Piilo J and Maniscalco S 2010 Frozen discord in non-Markovian depolarizing channels 

larXiv:1006.1805l 

[21] Mazzola L, Piilo J, and Maniscalco S 2010 Sudden transition between classical and quantum 
decoherence Phys. Rev. Lett. 104 200401 



15 



[22] Celeri L C, Landulfo AGS, Serra R M and Matsas G E A 2010 Sudden change in quantum and 

classical correlations and Unruh effect Phys. Rev. A 81 062130 
[23] Fanchini F F, Castelano L K and Caldeira A O 2010 Entanglement versus quantum discord in two 

coupled double quantum dots New J. Phys. 12 073009 
[24] Maziero J, Werlang T, Fanchini F F, Celeri L C and Serra R M 2010 System-reservior dynamics 

of quantum and classical correlations Phys. Rev. A 81 022116 
[25] Xu J S, Xu X Y, Li C F, Zhang C J, Zou X B and Guo G C 2010 Experimental investigation of 

classical and quantum correlations under decoherence Nat. Commun. 1 7 
[26] Xu J S, Li C F, Zhang C J, Xu X Y, Zhang Y S and Guo G C 2010 Experimental investigation of 

the non-Markovian dynamics of classical and quantum correlations Phys. Rev. A 82 042328 
[27] Soares-Pinto D O, Celeri L C, Auccaise R, Fanchini F F, deAzevedo E R, Maziero J, Bonagamba 

T J and Serra R M 2010 Nonclassical correlation in NMR quadrupolar systems Phys. Rev. A 

81 062118 

[28] Reina J H, Quiroga L and Johnson N F 2000 NMR-based nanostructure switch for quantum logic 
Phys. Rev. B 62 R2267 

[29] Khactskii A, Loss D and Glazman L 2003 Electron spin evolution induced by interaction with 

nuclei in a quantum dot Phys. Rev. B 67 195329 
[30] Imamoglu A, Knill E, Tian L and Zollcr P 2003 Optical pumping of quantum-dot nuclear spins 

Phys. Rev. Lett. 91 017402 
[31] Stepanenko D, Burkard G, Gicdkc G and Imamoglu A 2006 Enhancement of electron spin coherence 

by optical preparation of nuclear spins Phys. Rev. Lett. 96 136401 
[32] Lai C W, Maletinsky P, Badolato A and Imamoglu A 2006 Knight-ficld-enablcd nuclear spin 

polarization in single quantum dots Phys. Rev. Lett. 96 167403 
[33] Storcz M J and Wilhem F K 2003 Decoherence and gate performance of coupled solid-state qubits 

Phys. Rev. A 67 042319 

[34] Wang J, Batclaan H, Ponday J and Starace A F 2006 Entanglement evolution in the presence of 

decoherence J. Phys. B 39 4343 
[35] Scala M, Migliore R and Messina A 2008 Dissipation and entanglement dynamics for two 

interacting qubits coupled to independent reservoirs J. Phys. A 41 435304 
[36] Quiroga L, Rodriguez F J, Ramirez M E and Paris R 2007 Nonequilibrium thermal entanglement 

Phys. Rev. A 75 032308 

[37] Sinaysky I, Petruccione F and Bargarth D 2008 Dynamics of nonequilibrium thermal entanglement 
Phys. Rev. A 78 062301 

[38] Huang X L, Guo J L and Yi X X 2009 Nonequilibrium thermal entanglement in three-qubit XX 

model Phys. Rev. A 80 054301 
[39] Yan Y H, Wu C Q and Li B W 2009 Control of heat transport in quantum spin system Phys. Rev. 

B 79 014207 

[40] Sologubcnko A V, Lorcnz T, Mydosh J A, Rosch A, Shortslcevcs K C and Turnbull M M 2008 
Field-dependent thermal transport in the Haldane chain compound NENP Phys. Rev. Lett. 100 
137202 

[41] Sologubenko A V, Berggold K, Lorcnz T, Rosch A, Shimshoni E, Phillips M D and Turnbull M. M 
2007 Magnetothermal transport in the spin-1/2 chains of copper pyrazine dinitrate Phys. Rev. 
Lett. 98 107201 

[42] Lou P, Wu W C, Chang M C 2004 Quantum phase transition in spin-1/2 XX Heiscnberg chain 

with three-spin interaction Phys. Rev. B 70 064405 
[43] Krokhmalskii T, Derzhko O, Stolzc J and Vcrkholyak T 2008 Dynamic properties of the spin-1/2 

XY chain with three-site interactions Phys. Rev. B 77 174404 
[44] Yang M F 2005 Reexamination of entanglement and the quantum phase transition Phys. Rev. A 

71 030302(R) 

[45] Leggett A J, Chakravarty S, Dorsey A T, Fisher M P A, Garg A and Zwerger W 1987 Dynamics 
of the dissipative two-state system Rev. Mod. Phys. 59 1 



16 



[46] Breuer H P and Petruccione F 2002 The Theory of Open Quantum Systems (Oxford: Oxford 
University Press) 

[47] Goldman M 2001 Formal theory of spin-lattice relaxitiaon J. Magn. Reson. 149 160 

[48] Groisman B, Popescu S and Winter A 2005 Quantum, classical and total amount of correlations 

in a quantum state Phys. Rev. A 72 032317 
[49] Schumacher B and Westmoreland M D 2006 Quantum mutual information and the one-time pad 

Phys. Rev. A 74 042305 

[50] Wootters W K 1998 Entanglement of formation of an arbitrary state of two qubits Phys. Rev. Lett. 
80 2245 



